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i In the inflationary scenario all the matter constituting the universe was created from the 

Q\ , process of reheating after inflation. Recent development of the theory of reheating is briefly 

reviewed. The list of topics includes elementary (perturbative) theory of reheating; quantum 
field theory in a time- varying background; parametric resonance and explosive particle creation; 
non-thermal phase transitions from reheating; baryogenesis from reheating; residual oscillations 
of the scalar field, and other cosmological applications. 
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1. Introduction 



As a graduate student, I studied particle creation from strong gravitational and electric 
fields. I was impressed by a paper which Ya. Zel'dovich (1972) devoted to J. Wheeler's 
60th birthday volume. In particular, Zel'dovich quoted there the "Novikov paradox" : let 
an electron-positron pair be created from the strong electric field. The total momentum 
■ of the pair must be equal to zero, but the total energy is not zero. Query: how can 

this result be the Lorentz-invariant? Answer: the result indeed is the Lorentz-invariant, 
if the space-time interval between electron and positron is space-like. Hence, at the 
classical level, the creation of the pair is a causally independent event. The "Novikov 
• paradox" is a fair example of how Igor treats complex physical problems such as quan- 

pit turn creations of particles. Therefore taking occasion to contribute to Igor Novikov's 

60th birthday volume, I would like to consider another, one of the most spectacular, 
applications of the quantum theory of particle creation that emerges together with the 
Inflationary Scenario. Indeed, almost all matter constituting the Universe at the sub- 
sequent radiation-dominated stage was created from the reheating after inflation. 200z 
The term "reheating" here is an anachronism after the first inflationary models in which 
the Universe was hot before inflation and was reheated again after inflation. In modern 
versions of inflationary cosmology the pre-inflationary hot stage is no longer necessary 
(Linde 1990). It is assumed that the Universe initially expands quasi-exponentially in a 
vacuum-like state with a vanishing entropy and particle number density. At the stage 
of inflation, all energy was concentrated in a classical slowly moving inflaton field 4>. 
Soon after the end of inflation, when an observable universe was the size of a dime, the 
inflaton field began to oscillate near the minimum of its effective potential V(4>). An 
almost homogeneous inflaton field <j>(t) coherently oscillated with a very large amplitude 
of the order of the Planck mass cj> ~ M p . The interaction of the inflaton field with other 
elementary particles led to creation of many ultra-relativistic particles from the classical 
inflaton oscillations. Gradually, the inflaton field decayed and transferred all its energy 
non-adiabatically to the created particles. They interacted with each other and came to 
a state of thermal equilibrium at some temperature T r , which was called the reheating 
temperature. An enormous total entropy of the observable universe (or the total number 
of particles inside the horizon), S ~ 10 88 , was thus produced from the reheating after 
inflation. 

The reheating is an intermediate stage between the inflation and the radiation dom- 
inated Friedmann expansion. Therefore the reheating is associated with a number of 
issues of the Big Bang scenario: the large entropy problem, the baryogenesis problem, 
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the problem of relic monopolcs, the problem of topological defects, the primordial black 
holes problem, etc. 

The elementary theory of reheating based on the perturbation theory was developed 
right after the first models of inflation were suggested. Various aspects of the theory of 
reheating were further elaborated by many authors. The elementary theory of reheating 
found its place into the textbooks on the inflation (Linde 1990; Kolb & Turner 1990). 
Still, the general scenario of reheating in inflationary cosmology was absent, and a number 
of important problems remained unresolved. In particular, reheating in chaotic inflation 
theory, which covers many popular models of inflation, remained almost unexplored. 
Recently, the new effects in particle creation, arising beyond of the perturbation theory, 
were found, together with further constraints on the elementary theory itself. These 
effects significantly alter the reheating scenarios. In this contribution we briefly review 
the theory of the reheating after inflation. Wc will discuss the different approaches to 
the theory of reheating: from a quantitative, oversimplified point of view, and within 
rigorous qualitative formalism. This presentation is based on the ongoing collaboration 
with Andrei Linde and Alexei Starobinsky. 



2. Evolution of the Inflaton Field 

We consider a simple chaotic inflation. Among all the fields (</>, x, ip, At, hik, ■■■) which 
are present at these energy scales, the main contribution to gravity comes from the 
condensate of the scalar field </>, or from other condensates which effectively similar to 
that. In the fundamental Lagrangian C(((>, %, ip, Ai, hik, •••) f°r large energy density, we 
can retain for the cosmological applications only gravity and the dominant scalar field 

M 2 1 

C(R,4>) =--^R+-^-V{(l>). (2.1) 
1d7t 2 

The evolution of the FRW universe is described by the Einstein equation 

ff2 =5^(^ + H- (2 - 2) 

where H = a/a. The Klein-Gordon equation for (j>(t) is 

4> + 3H<p + V.^ = . (2.3) 



For a sufficiently large initial values of 4> > M p , a "drag" term 3Hip in (|2.3[) dominates 



over (j>; the potential term in (2.2) dominates over the kinetic term. This is the inflation- 
ary stage, where the universe expands (quasi-)exponentially, a(t) = aoexp(J dtH(tfj. 
However, with a decrease of the field (j> the "drag" term becomes less and less important, 
and inflation terminates at <f> < M p /2. After a short stage of the fast rolling down, the 
inflaton field rapidly oscillates around the minimum of V(<f>) with the initial amplitude 
<fio ~ O.lMp. Although this value is below the magnitude needed for inflation, it is still 
a very large figure. 

The character of the classical oscillations of the homogeneous scalar field depends on 
the shape of its potential V(<f>) around the minimum. We will consider two models: the 
quadratic potential V{<j>) = ^m^,(j) 2 and the potentia l V (4>) = j\(f> ' 



For the quadratic potential the solutions of eqs. ( |2.2| ) and ( |2.3| ) at the stage of oscil- 
lations are 

(j>{t) « cf> (t) ■ sin (m i) , 0o (t) = -=L ■ — . (2.4) 

V37T m^t 

The scalar factor averaged over several oscillations is a(t) w aot 2 / 3 . Oscillations of (j) in 
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this theory are sinusoidal, with the the amplitude 4>o(t) ~ 0.1M p a -3 / 2 decreasing as the 
universe expands. The energy of the field </> decreases in the same way as the density 
of nonrelativistic particles of mass m<f,: e^ = \<t> 2 + \m 2 h § 2 ~ a~ 3 . Hence the coherent 
oscillations of the homogeneous scalar field correspond to the the matter dominated 
effective equation of states with vanishing pressure. 

For the theory with the potential V(<f>) — jA0 it is more convenient to express the 
solutions of eqs. ( |2.2| ) and fl2.3| ) via the conformal time T] = J ^|y: 

4>(v) ~ Mv) ■ cn (^r), , 0o (77) = \J^ m p ■ ^ • ( 2 -5) 

The oscillations in this theory are not sinusoidal, but given by elliptic function 

cn( ^77, -j-), with the period of oscillations where a numerical constant c w 0.85 and 

the effective frequency of oscillations to = c^f\a<fio ~ 0.1\XM p . With a good accuracy 
one can write cj)(r)) ~ 4>oi r l) sin(c\/Aa0o'7). The scalar factor is a (t) oc r]{t) (x y/i. The 
amplitude of oscillations 4>o(t) ~ 0.1A/ p a _1 decreasing as the universe expands. The 
energy of the field </> decreases in the same way as the density of relativistic particles: 
e = \<t> 2 + 3^0 4 ~ a ~ 4 - The effective equation of state corresponds to the radiation 
dominated equation of states with pressure « |e^. 

3. Elementary Theory of Reheating 

The amplitude of the oscillations gradually decreases not only because of expansion of 
the universe, but also because of the energy transfer to particles created by the oscillating 
field. To describe the decay of the inflaton oscillations, we shall remember again the rest 
of the fundamental Lagrangian C(cj>, x, ip, Ai, hik, ■■■) which includes the other fields and 
their interaction with the inflaton field. The inflaton field cf> after inflation may decay into 
bosons x and fermions ip due to the interaction terms — ^g 2 <j) 2 X 2 an d —hipipcj), or into 
its own Bose quanta S(f) due to the self-interaction A(j) 2 5<fr 2 , or due to the gravitational 
interaction (like in the Starobinsky inflationary model). Here A, g and h are small 
coupling constants. In case of spontaneous symmetry breaking, the term — ^g 2 (j> 2 x 2 
gives rise to the term —g 2 a(f>x 2 - We will assume for simplicity that the bare masses of 
the fields x an d ip are very small, so that one can write m x (4>) « g<j), m^(<p) ~ \h<f>\. 

An elementary theory of reheating based on the perturbation theory was developed 
by Dolgov & Linde (1982) and by Abbot, Fahri & Wise (1982) for the new inflationary 
scenario. At the same time the theory of reheating was constructed in the Starobinsky 
(1982) model. Let us briefly recall the elementary theory of reheating. We consider for 
simplicity the classical field <fr(t) with the mass oscillating near the minimum of the 
quadratic potential. A homogeneous scalar field oscillating with frequency to — m$ can be 
interpreted as a collection of a number of 0-particles with zero momenta. The coherent 
wave of particles at rest with energy corresponds to particle density = e^/m^. 
Another way around, n$ oscillators of the same frequency m^, oscillating coherently 
with the same phase, can be viewed as a single homogeneous wave <fi(t). Based on that 
interpretation, the effects related to the particle production can be incorporated into 
the equation of motion ( |2.3[ ) for the inflaton field by means of the polarization operator 
(Linde 1990): 

4> + 3Hj>+ (ml + U{u)) = . (3.6) 
Here H(u>) is the flat space polarization operator for the field cj) with the four-momentum 
fci = (to, 0,0,0), to = rritj,. The real part of U-(to) gives only a small correction to m 2 , 
but when to > min(2m x ,2m^), the polarization operator H(to) acquires an imaginary 
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part Imll((jj). We will assume that 3> i? 2 ,m^ 3> Imll. The first condition is 
automatically satisfied after the end of inflation; the second usually is also true. The 



solution of (3.6) that generalizes the solution (2.4) and describes damped oscillations of 



the inflaton field is 



M p ( 1 



— — x exp --Tt sin (m^t) , (3.7) 
VoTrm<pt \ Z J 

where T is the total decay rate of ^-particles. Here we used a relation Im II = m$ • T, 
which follows from unitarity. Thus, eq. ( |3.7| ) implies that the amplitude of oscillations 
of the field <f> decreases as 4>q ( — ^rt) due to particle production 

which occurs during the decay of the inflaton field, as well as due to the expansion of the 
universe. 

For a phenomenological description of this effect one can just add an extra friction term 
T(f> to the classical equation of motion of the field </>, instead of adding the polarization 
operator (Kolb & Turner 1990): 

<j> + 3Htj> + T<j> + ml<f> = . (3.8) 



At the stage of oscillations the solution of this equation is similar to (3.7) 



One should note, that this equation, just as equation (3.7), is valid only under the 
conditions ^S> H, 3> V, and only at the stage of rapid oscillations of the field <j> 
near the minimum of V(<f>). This equation cannot be used to investigate of the stage of 
slow and fast rolling of the field </> during inflation. 

Suppose that T _1 is much less than the typical expansion time of the Universe 
the energy density of the inflaton field decreases exponentially within the time T : 
« 5 to2 j ( / , o ' exp(— Ti). This is exactly the result one would expect on the basis of 
the interpretation of the oscillating field <fi as a coherent wave consisting of decaying 
0-particles. 

The rate of decrease of the energy of oscillations coincides with the decay rate of <p- 
particles which is well-known from the perturbation theory. The rates of the processes 
4> ~ * XX an d 4> — > ipip (f° r m 4> ^ 2m x , 2m^) plus gravitational decay are given by 

Reheating completes at the moment t r when the rate of expansion of the universe given 
by the Hubble constant H = ^Jf^fz ~ t~ 1 becomes equal to the total decay rate V. The 
energy density of the universe at the moment t r — T~ 1 is 

3r 2 M 2 

<tr) - ■ (3.10) 

If the particles produced by the decaying inflaton field interact with each other strongly 
enough, then the thermodynamic equilibrium sets in quickly after the decay of the infla- 
ton field, and the matter acquires a temperature T r . The energy density of the Universe 
dominated by the ultrarclativistic particles in a state of thermal equilibrium is 

e(T r ) ~ ^T r 4 , (3.11) 
where a numerical factor g*(T r ) ~ 10 2 — 10 3 depends on the number of ultra-relativistic 



degrees of freedom. Comparing ( 3.10 ) and ( 3.11 ) we get the following result for the 
reheating temperature: 

T r ~ O.ly/T M p . (3.12) 
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Note that T r does not depend on the initial value of the field (j), it is completely determined 
by the parameters of the underlying elementary particle theory via the decay rate T. 

In order to get numerical estimates for the duration of reheating t r , for the decay rate of 
the inflaton field and for the reheating temperature T r , one should know the mass of the 
inflaton field and the coupling constants. We shall take into account several constraints 
on these in the framework of the perturbation theory (Kofman, Linde & Starobinsky 
1996b). The coupling constants of interaction of the inflaton field with matter cannot 
be too large, otherwise the radiative corrections alter the shape of the inflaton poten- 
tial. The necessary condition for the decay of inflaton oscillations is lo > 2m x ,2ro^,. 
Parameters of the inflaton potential are restricted from the constraints on amplitude of 
the cosmological fluctuations. All together it allows us to update a constraint on the 
elementary theory: the largest possible total decay rate in the perturbation theory is 
T < 10~ 20 M p . For instance, for the quadratic inflaton potential, it takes at least 10 14 
oscillations to transfer the energy of inflaton oscillations into the created particles. This 



is a very strong condition, which makes reheating very slow. From (3.12) one can obtain 



the general bound on the reheating temperature in the model of slow reheating: 

T r <10 9 GeV. (3.13) 

This is a very small temperature, at which the standard mechanism of baryogenesis in the 
GUTs cannot work; in this case we definitely need a theory of a low-temperature baryoge- 
nesis. At such a temperature no cosmologically interesting heavy strings, monopoles and 
textures can be produced, because the generation of the GUT topological defects require 
the high temperature phase transitions at Tqjjt ~ 10 16 GeV . For the massive inflaton 
field, in the slow reheating scenario the post-inflationary matter dominated equation of 
state of the inflaton oscillations lasts sufficiently long. It can enhance the gravitational 
instability of density fluctuations, which leads to the formation of the primordial black 
holes for some specific spectra of the initial density fluctuations. 



4. Quantum Field Theory in a Time-dependent Background 

Surprisingly, it was found recently that typically the transition between inflation and 
the hot Big Bang universe may be very different from what the elementary theory pre- 
dicts (Kofman, Linde & Starobinsky 1994; Shtanov, Traschen & Brandenberger 1995; 
Boyanovsky et al 1995a; Yoshimura 1995). The elementary theory of reheating might be 
still applicable if the inflaton field can decay into fermions only, with a small coupling 
constant h 2 <C m^/Mp. This theory also can provide a qualitatively correct description 
of particle decay at the last stages of reheating. However, the perturbation theory is 
inapplicable to the description of the first stages of reheating, which makes the whole 
process quite different. In what follows we will primarily consider the theory of the first 
stages of reheating. We will begin with the theory of a massive scalar field <f> decaying 
into particles x, then we consider the theory j4> 4 . 

For simplicity, we consider here the interaction — hg 2 F((f>)x 2 between the classical 
inflaton field ip and the quantum scalar field x with the Lagrangian: 

C(X) = \ (x,kX k ~ m\x 2 + tRx 2 ~ 9 2 FWx 2 ) ■ (4.14) 

Here £ is a coupling constant of interaction with the space-time curvature R. We also will 
distinguish two cases, the four-legs interaction g 2 F((f>)x 2 — 9 2 <t> 2 X 2 i an d the three-legs 
interaction g 2 F((f>)x 2 = 2g 2 o~4>X 2 which arises when symmetry is broken. The Heisenberg 
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representation of the quantum scalar field x 1S 

*(*• *) = (^372 / d " k (W*) + - ( 4 - 15 ) 

where and are the annihilation and creation Bose-operators. For the flat Fried- 
mann background with a scalar factor a(t) the temporal part of the eigenfunction with 
momentum k obeys the equation 

Xk + S^Xk+ + m 2 x - £R + g 2 F{4,)^Xk =0 . (4.16) 

Let us use the conformal time r\ and introduce the new function f k (rj) = a,(r\)Xk{j\) • Then 
instead of ( 4.16| ) we have 

/^ + n2/* = o, (4.17) 

where prime stands for and fi 2 = (m x a) 2 + k 2 + g 2 a 2 F((f>) + a 2 (l/6 — £)R. This equa- 
tion describes the oscillators with a variable frequency fl 2 , (t) due to the time-dependence 
of the background field (f>(rj) and 0,(77). At 77 — > -co we choose the positive-frequency 
solution /fc(»7) ~ e . We expect the quantum effect of the ^-particles creation and 
vacuum polarization as the inflaton field is varying during the slow- and fast-rolling down 
and oscillations after inflation. The problem is to find a general solution of the equation 



( 4.17 ) and regularize the formal expressions for the values of the energy-density (e x ) and 
the vacuum expectation (x 2 )- 



We adopt a physically transparent method to treat the eq. ( 1.17 ) for an arbitrary time 
dependence of the classical background field which was originally developed by Zeldovich 
and Starobinsky (1972) for the problem of the particle creation in varying strong grav- 
itational field. In terms of classical waves of the x _ne ld, quantum effects occur due to 
the departure from the initial positive-frequency solution. Therefore one may repre- 



sent solutions of eq.(4.17) as a product of its solutions in the adiabatic approximation, 



exp (±i J dig fife), multiplied by some functions 01(17) an d ftiv)'- 

f k (rj) = e -*fn kdv + hjri) e +iju k d v (4 lg) 

\/2f2fc \/2Cl k 
An additional condition on the functions a and (3 can be imposed by taking the derivative 



of the expression (4.18) as if a and (3 would be time-independent. Then the expression 
(4.18) is a solution of equation ( [1.17 ) if the functions a k ,f3 k satisfy the equations 



= J^-e+*/ n **> ft , fa = "L-M/n**, ak . (4.19) 

Normalization gives \a k \ 2 — \(3k\ 2 = lj and the initial conditions at rj — > —00 arc a k = 1, 
/3fe = 0. The coefficients a k (r]) and /3 k (rj) in our case coincide with the coefficients of the 
Bogoliubov transformation of the creation and annihilation operators, which diagonalizcs 
the Hamiltonian of the x _ne ld at each moment of time 77. 

The regularized vacuum expectation values for the x 2 > energy and particle number 
densities in terms of a k (rj) and (3 k (j]) are given correspondingly by 

(X 2 ) = m f°° dkk 2 ^-(\0 k \ 2 + Re(a k 0* k e- 2i I Q ^)) , (4.20) 

2lT 2 a 2 J Sl k \ ) reg 

1 / r 00 \ 1 



dkk 2 n k \/3 k \ 2 ) , (n x ) = -— dkk 2 \f3 k \ 2 . (4.21) 

/ ren Z7T d J n 



/0 ' reg 471 J Q 

The formal expressions for the vacuum expectation values should be renormalized. The 
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WKB expansion of the solution of eqs. (4.19) provides a natural scheme of regularization 



(Zel'dovich & Starobinsky 1972). Thus, the final result is expressed via the coefficients 
of the Bogoliubov transformation. 

The equation of motion of the inflaton field with the feedback of the quantum effects 

is 

+ 3i^ + y + i 3 2 (x 2 )F,0 = O , (4.22) 
where the vacuum average is given by eq.( 4.20| ). The quantum effects contribute to 



the effective mass of the inflaton field, e.g., for the massive scalar field rn 2 ^ 
^9 2 {x 2 )F,</>- The interactions with fermions —h'tpipcj) and the self-interaction —\<j) 2 5(j) 2 , 



where 6<f> < = <f> — <f>, result in extra terms h < ipip > and 3A<fi 2 (8<f) 2 ) in (4.22) 



With the self-consistent equations ( 4.20 ) and ( 4.22 ) based on the rigorous quantum 
field theory in the time varying background, one can investigate the role of the quantum 
effects on different stages of the evolution of the inflaton field. There are three main quan- 
tum effects in an external classical background: vacuum polarization, particle creation 
and generation of long-wavelength semi-classical fluctuations. Different effects dominate 
at different stages of the inflationary scenario. During inflation the long-wavelength fluc- 
tuations are generated for the minimally coupled light fields with £ — 0, but this effect 
is irrelevant to the reheating. The vacuum polarization and particle production could 
occur at the next stage of the fast rolling down of the inflaton field to the minimum of its 
potential V{<f>) at the end of inflation. However, these effects are negligible and do not 
alter the dynamics of the inflaton field. At the stage of oscillations of the inflaton field 
the only important effect is the creation of particles due to the non-adiabatic change of 
flk{v)- The particle production leads to the generation of (x 2 ), (e x ) and the damping of 
the inflaton oscillations. 

However, the term \g 2 {'X 2 )F ( j } in the equation (4.22) which describes the back re- 



action of the quantum effects on the evolution of the inflaton field is not reduced to 
the phenomenological "friction" term Tip, first suggested by Albrecht et al (1982). This 
oversimplified ansatz is widely used to investigate all the stages of the evolution of the 
inflaton field. Many authors used this method in order to understand whether reheating 
may slow down the rolling of the field <fi and, consequently, to support inflation. Remem- 
ber that one of the main reasons why inflation is possible is the presence of the friction 
term iHcj) in the equation of motion of the scalar field in an expanding Universe. The 
presence of a similar term due to reheating could support inflation even in the models 
where inflation would otherwise be impossible. As we seen, the damping of the scalar 
field during slow and fast rolling down cannot be described in such a simple way. In our 
study (Kofman, Linde & Starobinsky 1996b) we did not find any possibility to support 
inflation by reheating. As an exception, the friction term can correctly describe damping 
of oscillations in the case where the perturbation theory for interaction is applicable. 
It is instructive to return in the framework of the rigorous approach to the perturbation 



theory. Assuming small \(3k\ <C 1, from eqs. (4.19) one can obtain an iterative solution 



Pk^\ f H ^ exp {-2i I drf'Slktf'j) . ( 4 2:5 1 



A- 



Using expression for fife (77), we obtain that 0k{v) i s proportional to the integral 
9 2 J-oo drf e~ 2lkr > (a 2 F((j))y, where 4>(t) is given by the oscillating solutions ( [Q| ) or (2.5), 



and F{<j)) depends on the type of the interaction. This integral is evaluated by the method 
of stationary phase (e.g., Starobinsky 1982). To compare the results with these derived 
in the previous section, let us focus on the case of massive scalar field decaying via the 
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three-legs interaction g 2 <j<fix 2 - In this case the dominant contribution is given by the 
integration near r]k, where a{rjk) = This corresponds to the creation of a pair of 

^-particles with momentum k = ^a^/Jm^, from an inflaton with the four-momentum 
(m^, 0,0,0) at the moment r\ k . 

Substituting this estimation for (3k into (4.21), it is easy to calculate the decay rate 
of the inflaton fi eld. For the three-legs interactions guipx. or hfiipip it is indeed reduced 
to the formulas ( |3.9| ). For the four-legs interaction cfxf) — > XX, we have creation of a pair 
of x-particles with momentum k — a(rik)m^ from a pair of massive inflatons with the 
four-momentum (m^, 0,0,0) at the moment r]k- The decay rate of the massive inflaton 
field in this case is rapidly decreases with the expansion of the Universe as jr^(a 4e x) oc 
a~ 6 . Therefore the complete decay of the massive inflaton field in the theory with no 
spontaneous symmetry breaking or with no interactions with fermions is impossible! This 
observation gives rise to an unexpected possibility to consider the residual scalar field 
oscillations as a dark matter candidate. Indeed, if the inflaton oscillations correspond 
to the matter dominated equation of state, the inflaton field by itself, or other scalar 
fields with similar properties can be cold dark matter candidates, even if they strongly 
interact with each other. However, this possibility requires a fine tuning. More immediate 
application of our result is that it allows one to rule out the inflationary models which 
predict too small or too large value of the present radiation density parameter Q r . 

The fact that the number of particles |/3/c(?7)| 2 with the momentum k is predominantly 
generated at the certain moment r/k gives us a proper physical interpretation how the 
particles are created from the coherent oscillations of the homogeneous inflaton field: a 
pair of particles is created at the instant of the resonance k — uj ■ a(rj)/2 between k and 
the external frequency u) (which is equal for the massive inflaton). As the universe 
expands, that particular mode is redshifted from the resonance, and creation of particles 
with the momentum fc is terminated leaving n k — |/3fc(?7)| 2 <C 1. At an arbitrary moment 
r\ there is a corresponding momentum k resonating with the external frequency. 



5. Parametric Resonance: the Stage of Preheating 

A new interpretation that particles are created at the instant of resonance with the 
inflaton oscillations we have grounded in the last section gives us a hint for the general 
case when n k is not assumed to be small. 

Let us assume ther e is no symmetry breaking prior reheating. Therefore we will con- 
sider oscillations (2.4) of the massive scalar field <fi decaying into light particles \ due to 
the interaction — \g 2 4> 2 x 2 '-, an d oscillations (|2.5| ) in the t heory A</> 4 decaying due to the 



dHa^ Xk ) , ( e ,>...,„ n _ L x]l 



self- interaction — |a0 2 <5</> 2 . In the first case the equation (|4.16| ) for quantum fluctuations 
Xk{t) can be rewritten in the form 

J/2 v f , (j) f- +.9^o sin 2 (m^) + Aj (a :i ^ Xk ) = Q , (5.24) 

where <j>o stands for the amplitude of oscillations of the field </>, defined in the formula 
( [2.4| ), A = m 2 , + (9/2 — 3£)iJ 2 . As we shall see, the main contribution to x-particle 
production is given by excitations of the field x with k/a > m^, which is much greater 
than H at the stage of oscillations. Therefore, in the first approximation we may neglect 
the expansion of the Universe, taking a(t) as a constant and omitting the term A in 
( PI - 

In the case of theory with the j(j) 4 potential the role of quantum field plays the quan- 
tum fluctuations 8(f). In the equations ( 4.16| ) for the quantum fluctuations we have to 
substitute Xk — > S(pk, <? 2 — > 3A. Since in this theory the amplitude <fia(t) oc a~ 1 (t), it is 
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convenient to use the conformal time rj and equations for fluctuations in the form (4.17). 
For the oscillations (2.5) an exact equation for quantum fluctuations 5(f) can be reduced 
to the Lame equation. For simplicity here we will use an approximate equation 

d 2 (a5(j) k ) 



drj 2 



k 2 + 3Xa 2 (j)l sm 2 (cVXa(por]) (aSfa) = , (5.25) 



where c ~ 0.85, see (2.5). The crucial observation is that the equations (5.24) and (5.25), 
which describe the fluctuations in the models, can be reduced to the well-known Mathicu 
equation: 

d 2 y k 



dz 2 



(A(k) - 2qcos2z)y k = , (5.26) 



with the periodic effective frequency f2? = A(k) — 2gcos2z. For the first model (5.24), 
Uk — a 3 ^ 2 Xk, A(k) — -%s+2q, q — ^§r, z — cut, the frequency of the inflaton oscillations 
uj = m$. For the second model ( 5.25 ), y k = a- S(j) k , A w ^ + 2.08, and q w 1.04, z = lot], 
the frequency of inflaton oscillations uj = cy/~\a(f>o. 



An important property of solutions of the equation (5.26) is the existence of an ex- 
ponential instability y k oc expQi^z) within the set of resonance bands of frequencies 
Afc(") labeled by an integer index n. This instability corresponds to exponential growth 
of occupation numbers of quantum fluctuations nAt) oc exp(2fj,^ m^t) that may be in- 
terpreted as particle production. The simplest way to analyse this effect is to study the 
stability/instability chart of the Mathieu equation (MacLachlan 1961). 

We will distinguish different regimes of the resonance which correspond to the different 
regions of the parameters of the instability chart: 

Perturbation theory: q < 1 ; A — I 2 , I = 1, 2, ... ; 2it^k < H/w < 1 ; n k < 1 . If 
no expansion of the universe, fluctuations y k are slowly generated at the discrete modes k. 
The expansion quickly redshifts them from the resonance bands. The net effect is reduced 
to the creation of n k <C 1 particles with continuous spectrum, which we considered in 
the previous section. 

From this point of view the perturbative regime for the Starobinsky model was consid- 
ered by Starobinsky (1982) and for the chaotic inflation by Kofman, Linde & Starobinsky 
(1996b). 

Narrow resonance: q < 1 ; A ~ I 2 , I = 1, 2, ... ; 27r^ fc < H/w < 1 ; n k < O(l) . 
The universe expansion redshifts the momentum k/a from the narrow resonance band 
not too fast, so a nonsmall number of particles is generated at the resonance mode. This 
is an intermediate regime between the perturbation theory and the regime of the broad 
parametric resonance where expansion of the universe is a subdominant effect. 

The narrow resonance was mentioned by Dolgov & Kirilova (1990) for the new infla- 
tionary scenario. The importance of this regime for that model was first recognized by 
Traschen & Brandenberger (1990), but for several reasons their final results were not 
quite correct. A detailed theory of particle creation in the narrow resonance regime in an 
expanding universe was developed for the chaotic scenario by Kofman, Linde & Starobin- 
sky (1994; 1996b), see also Shtanov, Traschen & Brandenberger (1995) and Kaiser (1995). 
In order to investigate the explosive particle production from reheating, many authors 
considered the model of the narrow resonance without taking into account expansion of 
the universe because it made investigation simpler (Yoshimura 1995; Boyanovsky et al. 
1995a, b). However for some range of parameters many vital features of the reheating 
may disappear in this approximation. In particular, the effects in the model studied by 
Son (1996) disappear in an expanding universe. 

Broad resonance: q ~ 0(1) ; A > 2q ; 2irfi k ~ 0(1) ; n k 1 . Creation of 
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particles in the regime of a broad resonance is very different from that in the usually 
considered previous cases. In particular, it proceeds during a tiny part of each oscillation 
of the field <fi when 1 — cosz ~ q" 1 and the induced effective mass of the fluctuations 
y.k is less than m e tf. As a result, the number of the Bose particles grows exponentially 
fast. This regime occurs only if parameter q is nonsmall. A typical energy E of a particle 
produced at this stage is determined by equation A — 2q~ The line A = 2q divides 

the region of the broad resonance from the tachyonic regime where A < 2q ; fij. > 7r _1 . 
The line A — 2q also corresponds to the momentum k = 0. Near the line A — 2q typically 
/ifc ~ 0.175 in the instability bands, with the maximal value about 0.28. 

This regime is the most difficult to study. The investigation can be advanced with the 



following method: consider eq. (5.26) as the equation for the propagation of the wave 
y.k(z) backward in time. The evolution of yk{z) can be estimated for each scattering 
instance, then the net effect is the product of the successive scattering matrixes (Kofman, 
Linde & Starobinsky 1996b). The results for particle production based on this method 
were reported by Kofman, Linde & Starobinsky (1994). This stage also was studied by 
Boyanovsky et al. (1995a) and Fujisaki et al. (1995). 

We shall focus on the most interesting regime of the broad resonance. For our first 



model ( Ug) the necessary condition for it is -C gM p . A typical energy E of a 



particle produced at this stage is E ~ gm$M p . As we will see, the broad resonance 
(preheating) ends up within the short time t p h ~ 1 ln(m^ / g 5 M p ) . As a results, at the 
end of this stage the occupation number of created particles n& ~ exp{2pLkm c f ) t) with the 
energy k — E is extremely large: 

n E ~ \ > 1 . (5.27) 

For another model (5.25), the parameter q = 1.04 corresponds to the broad resonance 
regime. Looking at the instability chart, we see that the resonance occurs in the second 
band. The typical energy E of a created particles is E ~ Vkacfro. In the second band 
maximal value of the coefficient /j,k » 0.07. Note that the rigorous equation Lame for 
the fluctuations in the A</> 4 theory gives a twice smaller value of fik ~ 0.036. As long 
as the backreaction of created particles is small, expansion of the Universe does not 
shift fluctuations away from the resonance band, and the number of particles nk ~ 
exp(^at). The broad resonance in this model ends up within the time interval t p h ~ 
A _1 / 2 | In A|. At the end of this stage the occupation number is 

n E ~ \ > 1 . (5.28) 
A 

In addition to this models, the field 4> i n the theory A</> 4 may decay to x-particles due 
to the interaction —^g 2 <p 2 x 2 - This is the leading process for g 2 A. The parametric 
resonance is broad. The values of the parameter /i^ along the line A = 2q do not change 
monotonically, but typically for q ~ (? 2 /A ^> 1 they are 3 to 4 times greater than the 
parameter fik for the decay of the field <p into its own quanta. Therefore, the resonance 
in this theory is very efficient. 

The stage of broad parametric resonance which leads to the explosive particle creation 
is lasting typically for a few dozens of oscillations (for several oscillations in the last 
model). However, this does not mean that the process of reheating has been completed. 
Instead of created particles in the thermal equilibrium, by the end of this stage one has 
particles far away of equilibrium but with extremely large mean occupation numbers 
( |5.27 ), (5. 28). It is therefore convenient to divide the whole process of reheating into 



three different stages. At the first stage, which cannot be described by the elementary 
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theory of reheating, the classical coherently oscillating inflaton field <f> decays into bosons 
due to parametric resonance. In many models the resonance is broad, and the process 
occurs extremely rapidly (explosively). Because of the Pauli exclusion principle, there is 
no explosive creation of fermions. To distinguish this stage from the subsequent stages 
when particles decay or interact, we call it pre-heating. The second stage is the evolu- 
tion of the system after preheating. During this stage the particles produced from the 
preheating, can decay into other particles and self-interact. The last stage is the stage 
of thermalization, which can be described by standard methods. 



6. Back Reaction of Created Particles 

Creation of particles leads to the several effects which can change the dynamics of 
the system. In particular, it terminates the broad resonance particle production. First, 
the energy from the homogeneous field cf>(t) is transferred to the created particles. The 
amplitude of the classical oscillations cf>o is therefore decreasing faster than it woul d 
decrease due to the expansion of the universe only. Second, as it follows from ( 4.22 ), 
fluctuations contribute to the effective frequency of the inflaton oscillations: w^j = 
m + 9 2 (x 2 ) or w e// = c 2 X(f> 2 + 3A( 5<p 2 ). Third, fluctuations also can contribute to the 
effective mass of fluctuations themselves. As a result, the q parameters in the Mathicu 
equation are changing, in the different models correspondingly as q = g 2 ^/ '4w 2 j / and 
q = 3Aa 2 0§/4Wgy^. The A(k) parameters are changing as A(k) — k 2 /uj 2 j^a 2 + 2q, 
and A(k) = (k 2 + 3Aa 2 ( 5(f) 2 ))/uj 2 ^ + 2q. It moves the resonance momenta k out the 
their initial positions on the instability chart. The parameter q is decreasing, which 
dynamically shifts the parameters towards the narrower resonance. 

As a result, the creation of particles in the broad resonance regime is terminated. A 
simple criterion of the end of the preheating is the condition that the energy density of 
produced fluctuations is of the same order as the energy density the scalar field at the 
end of preheating t p h- For the first model (x 2 ) ~ 4%i f° r the second model ( S(f> 2 ) ~ 4>q. 
From this the timing of the preheating t p h given above was estimated. 

After copious particle creation around the resonance mode E slowed down, the infla- 
ton field still continues to transfer its energy into the energy of fluctuations. Son (1996) 
noticed that self-interaction of the S(f> particles in the Xcj) 4 theory might be important 
since they have very large occupation number density ue ~ 1/A at the resonance mode. 
Indeed, the usual estimation for the relaxation time r for the diluted gas of the 6(j> parti- 
cles due to the self-interaction has limited relevance when the occupation number density 
of particles is large. The scattering of created particles, however, does not eliminate the 
preheating regime. The particles created from the resonance, in the momentum space 
occupy only very narrow shells corresponding to successive resonance bands. Most of 
the particles after preheating are located in the first shell of the radius E and the width 
~ Q.1E (additionally their occupation number is sharply decreases towards the edges 
of the shell). Two particles from this shell can be scattered inside the shell. The rate of 
redistribution of particles within the shell is indeed much larger than r — 1 , but it does 
not terminate the resonance effect since particles remain inside the resonance shell. Two 
particle from the shell can be scattered outside of the resonance shell. However, this 
process is much slower than the preheating time t pr . There are interesting possibilities 
that some of the particles from the first resonance shell are scattered into other reso- 
nance shells. However, the whole process of rescattering and further creation of particles 
in the narrow resonance after preheating is rather complicated. Khlcbnikov & Tkachev 
(1996) performed first numerical simulation of the nonlinear evolution of the inflaton 
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field in A</> 4 theory with rescattering. Based on the result of Polarski & Starobinsky 
(1996) that the particles created with the large occupation number can be viewed as 
the superposition of the classical waves with random phases, they modelled preheating 
and scattering of <j) particles as the evolution of classical random gaussian field with the 
spectrum which has a sharp spike with the amplitude (5.28) at the mode E superim- 
posed with the homogeneous background. Numerical simulation confirms the features of 
preheating derived analytically. They also show the further growth of the fluctuations 
(64> 2 ) and reveal its spectrum after the preheating in this toy model. We shall note that 
in the more general case, when inflaton field can decay into other bosons, say, due to 
g 2 2 X 2 -interaction, the rescattering is less significant. We also do not find evidence of 
the strong Bose condensations at k = 0. 



7. Non-Thermal Phase Transitions from After Inflation 

As we have already seen, the reheating scenario depends on the type and strength of 
interactions of the involved fields. In this section we show that the reheating also strongly 
depends on the structure of the elementary particle theory. So far we considered theories 
with no symmetry breaking. In the case with symmetry breaking, in the beginning, 
when the amplitude of oscillations (f>o is much greater than ct, the theory of decay of the 
inflaton field is the same as in the case considered above. The most important part of 
pre-heating occurs at this stage. When the amplitude of the oscillations becomes smaller 
than m^/y/X and the field begins oscillating near the minimum of the effective potential 
at <fi = ct, particle production due to the narrow parametric resonance typically becomes 
weak. The main reason for this is related to the backreaction of particles created at the 
preceding stage of pre-heating on the rate of expansion of the universe and on the shape 
of the effective potential. However, importance of spontaneous symmetry breaking for 
the theory of reheating should not be underestimated, since it gives rise to the interaction 
term g 2 <j(f>x 2 or Xa(j)S(f> 2 which is linear in tfi. Such terms are necessary for a complete 
decay of the inflaton field in accordance with the perturbation theory. 

However, in the theories where preheating is possible one may expect many un- 
usual phenomena. One of the most interesting effects is the possibility of specific non- 
thermal post-inflationary phase transitions which occur after preheating (Kofman, Linde 
& Starobinsky 1996a, Tkachev 1996). These phase transitions in certain cases can be 
much more pronounced that the standard high temperature cosmological phase transi- 
tions ( Kirzhnits & Linde 1972, Linde 1990). They may lead to copious production of 
topological defects and even to a secondary stage of inflation after reheating. 

Let us first remember the theory of phase transitions in theories with spontaneous 
symmetry breaking. For simplicity let us consider the theory of scalar fields 4> the effective 
potential 

Vfax) = ^ -<**)* . (7.29) 

V{4>) has a minimum at (f> = cr, and a maximum at <f> = with the curvature = 
—m 2 = — Act 2 . This effective potential acquires corrections due to quantum (or thermal) 
fluctuations of the scalar fields AV = |A((<5</> 2 ))</> 2 . In the large temperature limit 
((54>) 2 ) = j^. The effective mass squared of the field (j> 

m\ eSf = -m 2 + 3A0 2 + 3A((<50) 2 ) (7.30) 
becomes positive and symmetry is restored (i.e. = becomes the stable equilibrium 
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point) for T > T c , where T 2 = ^j- to 2 . At this temperature the energy density of 

the gas of ultrarelativistic particles is given by p = g^(T c )^T* = 8m ^ffi^ • 

The theory of cosmological phase transitions is an important part of the theory of the 
evolution of the universe, and during the last twenty years it was investigated in a very 
detailed way. However, typically it was assumed that the phase transitions occur in the 
state of thermal equilibrium. Now we are going to show that similar phase transitions may 
occur even much more efficiently prior to thermalization, due to the anomalously large 
expectation values {(5(f>) 2 ) produced during the first stage of reheating after inflation. 

We will first consider the model ( [7.29 ) with the amplitude of spontaneous symmetry 
breaking a -C M p . As we seen in the previous sections, during preheating inflaton 
oscillates transfer most of its energy ~ AM* to its long- wave fluctuations {{8(f)) 2 ) ~ M 2 
in the regime of broad parametric resonance. 

The crucial observation is the following. If the initial energy density ~ \M p were 
instantaneously thermalized, the reheating temperature T r ~ A 1 / 4 Af p would be much 
greater than the typical particle energy after preheating E^ ~ \/AM p , and the magnitude 
of fluctuations ((S(f>) 2 ) ~ T 2 /12 ~ \f\~M 2 would be much smaller than the magnitude of 
non-thermalized fluctuations ((8(f)) 2 ) ~ M 2 . Thus after the first stage of reheating, the 
non-thermalized fluctuations of the scalar field <f> are much greater than the thermalized 
ones. Thermal fluctuations would lead to sym metry restoration in our model only for 
a < T r ~ X 1 / i M p . Meanwhile, according to eq. ( 7.3C| ), the non-thermalized fluctuations 
{(5(f)) 2 } ~ M 2 may lead to symmetry restoration in our model even if the symmetry 
breaking parameter a is as large as M p . Thus, the non-thermal symmetry restoration 
occurs even in those theories where the symmetry restoration due to high temperature 
effects would be impossible. Later on, ((6(f)) 2 ) oc a -2 , oc a -1 because of the expansion 
of the universe (as far as E$ ^> to). This leads to the phase transition with symmetry 
breaking at the moment t = t c ~ V~AM p m~ 2 when m^, e // = 0, {(8(f)) 2 ) — cr 2 /3, E x ~ to. 
Note that the homogeneous component (f>(t) at this moment is significantly less than 
\J ((6(f)) 2 ) due to its further decay after preheating. 

The mechanism of symmetry restoration described above is very general; in particular, 
it explains a surprising behavior of oscillations of the scalar field found numerically in 
the 0(iV)-symmetric model of Boyanovsky et al. (1995a). It is important that during 
the interval between preheating and the establishment of thermal equilibrium the uni- 
verse could experience a series of phase transitions which we did not anticipate before. 
For example, cosmic strings and textures, which could be an additional source for the 
formation of the large scale structure of the universe, should have a ~ 10 16 GeV. To 
produce them by thermal phase transitions in our model one should have the tempera- 
ture after reheating greater than 10 16 GeV, which is extremely hard to obtain (Kofman 
& Linde 1987). Meanwhile, as we see now, fluctuations produced at preheating may be 
quite sufficient to restore the symmetry. Then the topological defects are produced in 
the standard way when the symmetry breaks down again. In other words, production of 
superheavy topological defects can be easily compatible with inflation. 

On the other hand, the topological defect production can be quite dangerous. For 
example, the model (7.29) of a one-component real scalar field <f> has a discrete symmetry 
4> — ► —(f). As a result, after the phase transition induced by fluctuations ((8(f)) 2 ) the 
universe may become filled with domain walls separating phases (f> = +a and <f> = —a. 
This is expected to lead to a cosmological disaster. Investigation shows that in the 
theory ( [7.29 ) with a <C 10 16 Gev fluctuations destroy the coherent distribution of the 
background oscillations of <p and divide the universe in an equal number of domains 
with (f/ = ±<7, which leads to the domain wall problem. This mean that in consistent 
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inflationary models of the type of (7.29) one should have either no symmetry breaking 
or a > 10 16 Gev. 

Similar effects of the nonthermal phase transitions occur in the models where the 
symmetry breaking occurs for fields other than the inflaton field <p (Kofman, Linde & 
Starobinsky 1996a). The simplest model has an effective potential (e.g., Kofman & Linde 
1987; Linde 1994): 

Vfax) = \^ + f (x 2 - + ~<?W ■ (7.31) 

We will assume here that A <§; a < g 2 , so that at large c/> the curvature of the potential 
in the x-direction is much greater than in the (^-direction. In this case at large (f> the 
field x rapidly rolls toward x = 0. An interesting feature of such models is the symmetry 
restoration for the field x f° r 4* > 4>c — M/ g, and symmetry breaking when the inflaton 
field <f) becomes smaller than <fi c . As was emphasized by Kofman & Linde (1987), such 
phase transitions may lead to formation of topological defects without any need for high- 
temperature effects. 

Now we would like to point out some other specific features of such models. If the 
phase transition discussed above happens during inflation (i.e. if <j) c > M p in our model), 
then no new phase transitions occur in this model after reheating. However, for <j) c <C M p 
the situation is much more complicated. First of all, in this case the field <fi oscillates 
with the initial amplitude ~ M p (if M 4 < a\M p ). This means that each time when the 
absolute value of the field becomes smaller than (j) c , the phase transition with symmetry 
breaking occurs and topological defects are produced. Then the absolute value of the 
oscillating field <f> again becomes greater than <j> c , and symmetry restores again. However, 
this regime does not continue for a too long time. Within several oscillations, quantum 
fluctuations of the field x will be generated with the dispersion ((Sx) 2 ) ~ g~ x V\M 2 . For 
M 2 < g~ 1 V\aM 2 , these fluctuations will keep the symmetry restored. The symmetry 
breaking will be Anally completed when ((Sx) 2 ) will become small enough. Riotto and 
Tkachev (1996) noted that in the model ( 7.31 ) in the case of the strong self-interaction 
of x particles, A C j 2 « a, the presence of a huge occupation number of the created x 
particles ~ l/<? 2 would require the consideration beyond the one- loop approximation. 
It can be shown, however, that in this case the the preheating stops earlier when ~ 1/a 
(Kofman, Linde & Starobinsky 1996b). 

One may imagine even more complicated scenario when oscillations of the scalar field 
4> create large fluctuations of the field %, which in their turn interact with the scalar fields 
$ breaking symmetry in GUTs. Then we would have phase transitions in GUTs induced 
by the fluctuations of the field X- Note that in the models considered in this section the 
field x does not oscillate near x = prior to the phase transition, since such oscillations 
are damped out during the long stage of infl ation prior to the phase transition. Thus 
oscillations of the field x m the theory ( [7.31 ) definitely do not suppress the topological 
defect production. This means that no longer can the absence of primordial monopolcs be 
considered as an automatic consequence of inflation. To avoid the monopole production 
one should use the theories where quantum fluctuations produced during preheating are 
small or decoupled from the GUT sector. This condition imposes additional constraints 
on realistic inflationary models. 



8. Conclusions 



We have found that the decay of the inflaton field typically begins with the explosive 
production of particles during the stage of preheating in the regime of a broad parametric 
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resonance. During the second stage, the inflaton field decays further in the regime of 
the narrow resonance, and particles produced at this stage decay into other particles and 
self-interacting. The third stage is the thermalization. The two last stages require much 
more time that the stage of preheating. Typically coupling constants of interaction of 
the inflaton field with matter are extremely small, whereas coupling constants involved in 
the decay of other particles can be much greater. As a result, the reheating temperature 
T r given by (3.12) will not necessarily be defined by the perturbative decay rate of the 
inflaton oscillations (3.E). T r can be much higher than the typical temperature T r < 10 9 
GeV, which could be obtained neglecting the stage of parametric resonance. 

A specific feature of the preheating is that bosons produced at this stage are far away 
from thermal equilibrium and typically have enormously large occupation numbers. This 
may have very interesting applications to cosmology. There is a new class of phase 
transitions that may occur at the intermediate stage between the end of inflation and 
the establishment of thermal equilibrium. These phase transitions take place even in the 
theories where the scale of spontaneous symmetry breaking is comparable to M p and 
where the reheating temperature is very small. Therefore, phase transitions of the new 
type may have dramatic consequences for inflationary models and the theory of physical 
processes in the very early universe. 

The preheating may help with the baryogcncsis problem (Kofman, Linde & Starobinsky 
1996a; Yoshimura 1996). In the models of GUT baryogenesis, it was assumed that 
GUT symmetry was restored by high-temperature effects, since otherwise the density 
of X, Y, and superheavy Higgs bosons would be very small. This condition is hardly 
compatible with inflation. It was also required that the products of decay of these 
particles should stay out of thermal equilibrium, which is a very restrictive condition. 
In our case, the superheavy particles responsible for baryogenesis can be abundantly 
produced by parametric resonance, and they as well as the products of their decay will 
be out of thermal equilibrium until the end of reheating. 

Another consequence of the resonance effects is a fast change of the equation of state 
from a vacuum-like one to the equation of state of relativistic matter. This leads to 
suppression of primordial black holes formation that could be produced after inflation. 

The preheating stage may appear not only in the models where the reheating occurs 
due to the decay of the homogeneous scalar field. Kolb & Riotto (1996) suggested 
that preheating may also occur in first order inflation. It would be very interesting to 
investigate this problem, as well as the many other exciting possibilities which the new 
theory of reheating may offer us. 
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